第 41 卷 第 4 期 天 文 学 进展 Vol. 41, No.4 
2023 年 12 月 PROGRESS IN ASTRONOMY Dec., 2023 


doi: 10.3969 /j.issn. 1000-8349.2023.04.10 
基于 卡尔 曼 滤 波 与 Camshift 的 低 轨 空间 
变 光 目标 跟 踩 方法 


ERR, Preemie, ERA, EIR, oe 胜 5 


(1. 昆明 理工 大 学 信息 工程 与 自动 化 学 院 ， 昆 明 650031; 2. 中 国 科学 院 紫金 山 天 文 台 空间 目标 与 碎片 观测 
重点 实验 室 ， 南 京 210023; 3. 琼 台 师 范 学 院 信息 科学 技术 学 院 ， 海 口 571127; 4. 楚雄 师范 学 院 数学 与 计 
算 机 科学 学 院 ， 楚 雄 675099; 5. 华中 科技 大 学 人 工 智能 与 自动 化 学 院 ， 武 汉 430074) 


摘要 : 针对 低 轨 空间 变 光 目 标 星 像 形态 变化 频繁 且 运 动 速度 快 ， 跟 踪 目 标 时 很 难 准确 捕 绪 目标 
位 置信 息 这 一 难点 ， 提 出 卡尔 曼 滤 波 与 Camshift 相 结 合 的 目标 跟踪 方法 。 根 据 目 标的 难 捕获 4 
性 ， 引 入 卡尔 曼 滤 波 外 推 对 目标 进行 位 置 预测 ， 用 改进 的 Camshift 跟踪 目标 ， 即 基于 灰 度 图 
单 通道 的 白色 提取 的 目标 跟踪 方法 ， 实 现 低 轨 变 光 目 标 稳健 跟踪 和 提高 目标 捕获 率 ; 当 目标 与 
星 瞬 间 相互 遮挡 时 ， 用 改进 的 遮挡 目标 预 判 方法 进行 位 置 预测 ， 同 时 用 卡尔 曼 滤 波 的 预测 位 置 
 Camshift 计算 出 的 目标 位 置 ， 作 为 观测 位 置 去 更 新 卡尔 曼 滤波 ， 实 现 遮 挡 目 标 稳健 跟踪 和 
高 目标 捕获 率 。 实 验 结果 表明 : 跟踪 目标 的 可 变 波 门 根据 目标 大 小 自 适应 调整 ， 不 仅 可 以 实现 
高 轨 目 标 稳定 跟踪 ， 而 且 对 低 轨 变 光 目标 具有 稳健 跟踪 效果 。 当 目标 与 恒星 瞬间 相互 遮挡 时 
稳健 跟踪 目标 ， 提 高 跟踪 数据 的 有 用 率 。 该 方法 程序 运行 速率 快 ， 灵 敏 度 高 ， 适 用 性 强 ， 实 时 
好 ， 具 有 很 好 的 使 用 价值 和 广泛 应 用 前 景 。 

X 键 id: 卡尔 曼 滤波 ，Camshift， 变 光 目 标 ， 跟 踪 

中 图 分 类 号 : P114.4 文献 标识 码 : A 


q 


加 | N 


ma 
E Rot Hos 


1 引 F 


空间 目标 跟踪 是 获取 目标 轨道 信息 的 关键 手段 之 一 ， 跟 踪 的 鲁 棒 性 直接 关系 到 空间 目 
标 自 动 识别 、 定 位 和 编目 及 望远镜 运行 效率 。 随 着 国际 运动 目标 研究 技术 不 断 提高 ， 空 间 目 
标 跟 踪 方法 由 手动 跟踪 逐步 向 自动 跟踪 转变 ， 目 标 跟 踪 数 量 由 单 目 标 跟 踪 向 多 目标 跟踪 方 
向 发 展 ”。 目 前 跟踪 方式 正在 由 地 基 跟 踪 向 天 基 跟 踪 逐 步 发 展 ， 为 空间 态势 感知 开辟 更 广 
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阔 发 展 方向 和 应 用 情景 。 

尽管 空间 目标 跟踪 方法 不 断 改 善 ， 但 稳定 跟踪 低 轨 空间 变 光 目标 仍 是 一 个 挑战 性 难题 。 
首先 ， 低 轨 空 间 变 光 目标 运动 速度 较 快 ， 导 致 跟踪 波 门 跟 不 上 目标 运动 速度 ， 容 易 跟 丢 目 
标 ; 其 次 ， 目 标 与 恒星 存在 瞬间 相互 遮挡 ， 有 瞬间 遮挡 是 指 被 跟踪 目标 前 一 帧 或 后 一 帧 与 恒星 
没有 相互 遮挡 ， 而 当前 帧 与 恒星 存在 遮挡 ， 当 前 帧 目标 与 恒星 星 像 重 倒 变 成 一 颗 更 大 的 不 规 
则 目标 ， 致 使 目标 跟踪 误差 过 大 ， 甚 至 无 法 捕获 目标 位 置 ， 再 次 ， 固 定 距 ( 波 门 ) 方法 是 使 
最 广泛 的 计算 空间 目标 质心 方法 ， 但 变 光 目标 星 像 形态 变化 频繁 ， 目 标 圆 整 度 差 ， 准 确 计 
算 波 门 质心 较 困 难 ， 传 统 的 固定 波 门 跟踪 方法 不 适应 星 像 变 化 ， 且 捕获 到 较 少 的 位 置信 息 ， 
降低 了 目标 捕获 率 ， 问 题 有 待 解决 。 如 图 Ata, Bla), b) 和 c) 的 红色 方 框 为 跟踪 波 门 ， 
近似 白色 圆 点 为 空间 目标 ， 其 余 白 色 长 条 状 目标 为 恒星 。a) 图 中 波 门 跟 不 上 目标 运动 速度 
导致 目标 跟 丢 ; b) 图 中 目标 与 恒星 相互 遮挡 ， 没 有 捕获 位 置信 息 或 者 得 到 误差 大 的 位 置信 
息 ; c) 图 中 目标 比 波 门 大 ， 不 易 捕 获 目标 位 置信 息 。 


b) 
1 ， 低 轨 变 光 目 标 跟 踪 存 在 的 关键 问题 


目前 国内 外 主要 使 用 具有 普遍 适用 性 的 卡尔 曼 滤 波 跟 踪 空 间 目标 ， 科 学 家 们 也 基于 卡 
尔 曼 滤 波 、Camshift 和 深度 学 习 相关 算法 等 提出 一 些 空间 目标 跟踪 方法 。Chai 等 人 “及 任 
楠 “提出 用 卡尔 曼 滤波 和 深度 学 习 相关 的 空间 目标 跟踪 方法 来 减 小 跟踪 误差 。 平 一 蜀 " 提出 
空间 目标 与 恒星 在 相互 遮挡 情 时 的 星 像 分 离 方 法 。 王 均 和 许 占 伟 “、Li SEA BR 
出 的 方法 重 在 增强 跟踪 稳健 性 和 适用 性 。 许 占 伟 和 王 坎 ” 提 出 基于 轮廓 特性 的 目标 跟踪 ， 方 
法 侧重 于 解决 低 轨 空间 变 光 目 标 跟 踪 ， 但 算法 相对 复杂 且 未 提 及 变 光 目标 与 恒星 互相 遮挡 
的 处 理 方法 。 总 之 ， 上 述 空间 目标 跟踪 方法 虽然 具有 一 定 适用 性 ， 但 没有 一 个 方法 可 以 同时 
解决 低 轨 空间 变 光 目标 跟踪 面临 的 上 述 三 个 难题 。 

鉴于 卡尔 曼 滤 波 具 有 最 佳 位 置 预测 的 优点 ，Camshift 具有 自 适应 调整 目标 窗口 的 中 心 
位 置 与 大 小 的 优点 ， 且 程序 运行 速度 快 、 实 时 性 好 ， 同 时 在 空间 目标 跟踪 时 ， 常 将 窗口 质心 
作为 目标 质心 。 为 此 将 两 者 结合 ， 以 解决 低 轨 空间 变 光 目标 跟踪 面临 的 三 个 难题 。 尽 管 卡尔 
曼 滤 波 与 Camshift 的 结合 曾 应 用 于 汽车 、 轮 船 、 行 人 等 运动 目标 跟踪 ， 但 目前 国内 外 还 没 
有 将 其 应 用 于 空间 变 光 目标 跟踪 的 报道 ， 本 文 首次 应 用 于 低 轨 空间 变 光 目标 跟踪 。 另 外 ， 不 
同 于 传统 的 卡尔 曼 与 Camshift 的 结合 ， 本 文 针 对 低 轨 变 光 目 标 跟踪 问题 提出 基于 卡尔 曼 滤 
波 与 改进 的 Camshift 结合 的 跟踪 方法 ， 同 时 提出 目标 瞬间 遮挡 的 位 置 预 判 和 跟踪 方法 ， 获 
得 稳健 跟踪 效果 。 
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2 卡尔 曼 滤 波 与 Camshift 的 混合 算法 


如 图 2 所 示 ， 卡 尔 曼 与 Camshift 结合 的 混合 算法 主要 分 三 个 阶段 : 首先 用 背景 减 除 
法 检测 出 空间 目标 ; 其 次 用 卡尔 曼 滤波 对 空间 目标 进行 位 置 预 测 ; 然后 根据 预测 结果 用 
Camshift 跟踪 目标 。 当 用 改进 Camshift 跟踪 目标 时 ， 分 低 轨 变 光 目标 、 中 高 轨 目 标 和 目标 
遮挡 三 种 情况 。 针 对 变 光 目标 和 中 高 轨 目 标 ， 在 目标 与 恒星 互 不 遮挡 时 ， 用 改进 Camshift 
跟踪 目标 和 捕获 目标 位 置 ， 针 对 目标 与 恒星 互相 遮挡 ， 用 遮挡 目标 跟踪 方法 来 跟踪 目标 
和 捕获 目标 位 置 ， 最 后 更 新 卡尔 曼 滤 波 器 ， 预 测 下 一 帧 目标 位 置 ， 然 后 返回 到 卡尔 曼 与 


Camshift 结合 的 跟踪 目标 步 又， 如 此 反复 进行 。 


NI 


| i T 
改进 的 Camshift 跟踪 目标 遮挡 目标 跟踪 方式 
| 


2 ”卡尔 曼 与 Camshift 结合 跟踪 目标 


2.1 ”卡尔 曼 滤 波 预测 目标 位 置 

令 状 态 向 量 Xp = [zk, Yk, Vek, vyl” 的 各 个 值 分 别 代表 空间 目标 z A y 方向 的 位 置 和 
速度 ， 观 测 向 量 Ze = [zx, yx)? 代表 空间 目标 的 观测 位 置 ， 观 测 位 置 由 Camshift 跟踪 目标 
时 获得 。 其 状态 方程 和 输入 的 观测 方程 分 别 为 : 


Xk = ÅkXk-ı + BkWk , (1) 
Zr = HkXk +t Ve, (2) 


HR, X, 和 Xp 分 别 为 和 一 1 时刻 的 状态 向 量 ，Ak 为 状态 转移 和 矩阵， 万 为 输入 
JERE, A, 为 观测 矩阵 。Wi 和 Wi 分 别 为 状态 噪声 和 观测 噪声 矩阵 ， 对 应 的 方差 矩阵 分 别 
AQ 和 RR。 其 推导 公式 如 下 : 


4 期 


HP, Xk AKA TUE, Pl SWIA, Ke A 


值 ， Pk 为 Xk 协 方差 矩阵 。 


(3) 和 (4) 式 的 猜测 
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Xj = AkXpr-1 + BaeWh_1 , (3) 
P = AgPp1A'+Q , (4) 
Ky =P, HE (HP; He +R) , (5) 
X, =X, + Ki(Ze— HX) ， (6) 
P, = (I — Kh Hp) P; , (7) 


FE 尔 曼 增益 矩阵 ，Xh AX 的 修正 


立 置 不 是 最 佳 的 预测 位 置 ， 因 为 没有 考虑 当 


前 时 刻 Camshift 跟踪 目标 时 获得 的 观测 位 置 ， 通 过 (5)、(6) 和 (7) 式 对 当前 的 观测 位 置 进 
行 更 新 修正 后 输出 最 佳 预 测 位 置 ”“。 
2.2 Camshift 跟踪 目标 
Camshift 根据 卡尔 曼 滤 波 获得 最 佳 预测 位 置 对 空间 目标 跟踪 。 在 Camshift 算法 中 ， 设 
T(x, Yk) 是 坐标 (£k, Yk) 像素 点 的 像素 值 ， 搜 索 窗 fà H 的 零 阶 距 Moo> 一 阶 距 为 Myo. Moi» 
则 搜索 窗口 质心 tko Yk 和 更 新 搜索 窗口 S 的 大 小 分 别 为 : 
_ Mio 
Tke = Moo > (8) 
Moi 
Yke Moo , (9) 
E Moo 
5 =) 3 (10) 
民 据 Moo 调整 搜索 窗口 大 小 ， 并 将 搜索 窗口 中 心 逐渐 移动 到 质心 当 移 动 距离 大 于 预 
设 固 定 闵 值 ， 重 新 计算 调整 后 的 窗口 质心 ， 直到 窗口 中 心 与 质心 之 间 的 移动 距离 小 于 预 设 固 
定 阔 值 ， 或 者 循环 运算 次 数 达到 最 大 值 ， 则 认为 收敛 ， 进 入 下 一 帧 的 目标 搜索 。 


2.3 ”改进 的 Camshift 


Camshift 运行 效率 高 ， 对 CCD 图 
目标 为 白色 。 传 统 的 Camshift 


BEI 
黑色 ， 


行 目标 跟踪 ， 色 


时 的 HSV # 
ZELMER 


0 颜色 对 比 度 低 等 性 


pa 


Eg 


背 [三 
月 是 


色 度 万 的 值 为 0" 一 360"*， 饱 和 
色 度 万 需要 对 红 绿 蓝 等 多 个 颜色 
存在 噪点 、 复 杂 背 景 


k 3 个 通道 提取 颜 
SULT AA BR 


:下 的 跟踪 效果 较 好 ， 其 中 CCD ae 
法 将 RGB 颜色 转换 为 HSV 颜色 并 根据 颜色 信息 
E 5 的 值 为 0~1， 亮 度 V 的 值 为 0 一 1， 


色 


jo 


显然 颜色 提取 是 一 个 复杂 过 程 ， 当 
踪 错 误 ， 为 此 对 Camshift 跟踪 目标 


Bl 


IZ 


跟踪 目 


GRAY 


提取 


对 的 最 佳 


标 时 


色 模 型 进行 改进 。 


色 里 的 


Ñ H 


0. S 0. V 


1。 但 空 


由 于 空间 目标 图 像 是 一 种 
颍 色 提 取 改 为 基于 单 通道 gray 图 
颜色 提取 的 复杂 度 。 

图 像 在 RGB Pil 


像 的 白 1 


实现 原理 如 下 。 
AEN R= G= B, GRAY 


色 提 取 并 跟踪 


通道 的 灰 度 图 像 ， 于 是 将 HSV 
标 ， 这 样 可 以 大 大 简化 


色 


C4 


色 在 HSV ği 


图 像 的 


间 目 标的 


白色 不 一 定 是 100% 的 纯 


色 ， 为 此 
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设置 饱和 度 S 的 最 大 值 为 P， 亮 度 V 的 最 小 值 为 Q, 


此 处 忆 的 值 为 0.25， 


H=0, 
0<S<P 


JVI. 


Q 的 值 为 0.75 比较 适宜 。 


2.4 ”改进 的 遮挡 目标 预 判 和 跟踪 方法 


Camshift 结合 的 方法 具 


则 颜色 提取 算法 改 为 如 下 : 


目标 与 恒星 瞬间 氮 挡 在 低 轨 变 光 目 标 跟踪 中 是 普 裔 存在 的 现象 ”将 卡尔 曼 滤 波 与 


有 根据 目标 大 小 


但 目标 与 恒星 瞬 


时 互相 i 


pe 


ASS 


HERREY 


2.4.1 遮挡 目标 预 判 
(1) 
连 近 恒星 预 判 非 
离 预 判 ， 


与 恒星 质心 之 间 的 距离 巴 


寻找 跟踪 目标 周围 的 连 近 恒星 


常 重要 ， 如 果 对 整 帧 图 


判 ， 则 无 法 准确 


内 目标 与 4 


面积 9 为 : 


吾 星 之 间 的 距离 预 判 ， 


= |n(max(z;) 


H Pr- E Æ HA 


时 ， 由 于 恒星 的 干扰 无 法 准 


自 适应 调整 目标 波 门 和 准确 
时 ， 将 该 处 目标 与 恒星 当 作 
从 位 置 与 实际 位 置 之 间 的 误差 进 
目标 的 相应 解决 办 法 。 


像 进行 大 范围 
判断 目标 与 恒星 的 遮挡 
被 跟踪 目标 质心 为 矩 的 质心 


E 标 和 纵 坐 标的 扩大 倍数 ， 


— min(z:))| - In(max(y) 


wait § 


KAR YE YS 3 


需要 1 
一 


行 位 置 
帧 恒星 的 实际 位 置 来 预测 当前 帧 的 可 


预测 ， 


的 目标 质心 与 恒星 质心 之 间 的 距 


预测 目标 位 置 的 优点 。 
一 颗 更 大 目标 跟踪 ， 使 得 该 帧 被 跟踪 
步 增 大 ， 失 去 跟踪 价值 ， 为 此 提出 针对 低 轨 变 光 


必然 增加 计算 量 和 耗 时 ， 降 低 程序 运行 效率 ;如 果 只 对 目标 周围 小 范 


大 的 目标 质心 


问题 。 


FRERES H 


标 周 


— min(yi)) 


此 处 采用 特定 矩 面 积 
， 被 跟踪 目标 波 门 的 横 坐 标 为 
围 连 近 恒星 的 矩 


(14) 


出 目 


能 位 置 ， 即 月 


标的 实际 质心 位 置 ， 这 时 
用 上 一 帧 被 跟踪 目标 的 实际 位 置 来 预测 当前 帧 的 可 能 位 置 


~ 


预测 位 置 代 蔡 实际 位 置 进行 目标 


与 恒星 质心 之 


By 


通常 采用 搜索 窗 


波 预 测 出 目标 当前 


间 的 位 置 预 
(3 ) 目标 与 恒星 的 预 党 


贞 与 连 近 恒星 之 间 的 距离 。 


ll > 


I 位 置 之 间 的 距离 计算 
口 的 质心 作为 空间 目标 


当 有 多 帧 目标 追 挡 时 需要 进行 多 


的 质心 和 连 近 恒星 的 质心 ， 
设 当前 帧 被 跟踪 目标 的 预测 质心 为 


贞 预 测 。 


星 是 否 存在 互相 遮挡 ， 将 目标 与 恒星 之 间 的 预测 距离 设置 阔 值 


连 近 恒星 的 预测 质心 为 (xi, we)j， 则 目标 与 恒星 之 间 的 预测 距离 万 为 ; 
D= V (Lke 一 Lie) 2 十 (Yke = Yie)? 
(4) DAE Be 
为 更 好 地 判断 目标 与 恒 
F, BU Zao F: 
i D<F 
K= 
0, DSF 


然后 通过 卡 


iE. A> = yy 
ROR Sve 


(Bras Uke)» 


(15) 


(16) 
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当 恒 星 与 目标 质心 之 间 的 预测 距离 D KAFEER FIN, RK 值 为 0， 表 示 不 存在 目 
标 遮 挡 ， 当 目标 与 恒星 之 间 的 预测 距离 D 小 于 给 定 闵 值 时 ， 取 KK 值 为 1， 表 示 存 在 目标 
遮挡 。 
2.4.2 遮挡 目标 跟踪 

当 预 判 存在 遮挡 目标 时 ， 用 卡尔 曼 滤 波 的 预测 位 置 代替 Camshift 算法 计算 出 的 目标 
位 置 作为 观测 值 去 更 新 卡尔 曼 滤波 。 遮 挡 目 标 跟 踊 结 束 时 ， 跟 踪 方 式 又 切换 到 卡尔 曼 与 
Camshift 结合 的 跟踪 算法 。 


Kir 


7IT 


3 x 验 


实验 以 中 国 科 学 院 紫 金山 天 文 台 姚 安 观测 站 500 mm 口径 光学 望远镜 在 2021 年 2 月 
4 日 拍摄 的 多 帧 运动 目标 图 像 为 素材 。 基 于 戴尔 笔记 本 电脑 WINDOWS10 系统 上 使 用 
pycharm2022 软件 ， 基 于 python3.9 的 opencv3.4.17.63 程序 。 固 定 波 门 跟踪 用 空间 目标 观 
测 中 常用 的 卡尔 曼 滤波 ， 自 适应 的 可 变 波 门 跟踪 将 卡尔 曼 滤 波 与 Camshift 相 结合 。 跟 踪 精 
度 计算 通过 计算 跟踪 误差 的 RMS， 而 RMS 是 跟踪 目标 时 获得 的 估计 值 与 预测 值 之 差 。 
3.1 Camshift 跟踪 目标 
3a) Alb) 两 帧 图 中 的 红 、 绿 、 蓝 三 个 波 门 分 别 代表 测量 波 门 、 预 测 波 门 、 状 态 波 门 。 
a) 用 Camshift 跟踪 目标 ， 一 般 跟 踪 稳定 ， 偶 尔 存 在 目标 跟踪 异常 的 情况 发 生 。b) 用 改进 
Camshift 跟踪 目标 ， 方 法 能 稳健 提取 目标 颜色 信息 和 跟踪 目标 ， 三 个 波 门 收 敛 于 目标 质心 ， 
在 相同 运行 环境 下 与 Camshift 跟踪 相 比 ， 程 序 的 运行 效率 提高 了 约 1.17 倍 。 


X 


a) b) 
3 Camshift 跟踪 目标 


表 1 对 多 圈 被 跟踪 目标 的 异常 情况 进行 统计 ，Camshift 跟踪 目标 时 的 平均 异常 率 为 
4.6%, zt Camshift 跟踪 目标 时 无 异常 情况 发 生 。 


表 1 跟踪 目标 异常 对 比 


方法 采集 时 间 ”目标 ”平均 每 圈 ”平均 每 圈 跟 ”异常 率 
/min 圈 数 | 采集 帧 数 ” 踪 异常 帧 数 /% 
Camshift 3 17 65 3 4.6 


改进 Camshift 3 17 65 0 0 
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3.2 ”中 高 轨 目 标 跟踪 


Al 4 
图 4a). b) c) 用 卡尔 曼 滤 波 的 固定 波 门 跟踪 中 高 轨 目 标 ，a) 中 目标 大 于 波 门 ，b) 和 
c) 中 的 目标 小 于 波 门 ， 该 方法 能 跟踪 目标 ， 但 波 门 不 具有 收敛 性 。d)、e)、f) 三 帧 图 改 用 卡 
尔 曼 滤波 与 Camshift 结合 的 自 适 应 可 变 波 门 跟踪 中 高 轨 目 标 ， 波 门 收敛 于 目 


标 质 心 。 


d) e) 


图 4 


固定 波 门 与 可 变 波 门 跟踪 中 高 轨 目 标 
#2 的 卡尔 曼 与 Camshift 结合 的 中 高 轨 空 间 
获 率 ， 目 标 捕 


目标 捕获 率 略 高 于 卡尔 曼 跟 踪 时 的 目标 捕 
捕获 率 平 均 提 高 6%。 
表 2 中 高 轨 目 标 捕获 率 对 比 
方法 采集 时 间 标 ”平均 每 圈 ” 波 门 大 小 是 否 “平均 每 圈 位 “捕获 率 
min Až ”采集 帧 数 。” 自 适 应 变化 。 置 捕获 帧 数 /% 
卡尔 曼 3 15 67 T 58 86.5 
卡尔 曼 与 Camshift HAA 3 15 67 是 62 92.5 
图 5 的 卡尔 曼 与 Camshift 结合 的 跟踪 精度 优 于 卡尔 曼 的 跟踪 精度 ， 跟 踪 精 度 平 均 提 高 
了 1.21 倍 ， 获 得 较 好 的 跟踪 效果 


目标 编号 
, ERA KK SS AGH Camshift 4 


123 4 5 6 7 8 9 10 1112 13 1415 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 
E: HERR ES, BE 


图 5 


中 高 轨 目 标 精 度 统 计 
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3.3 
3.3.1 


下 角 ， 


卡尔 曼 与 Camshif 结合 ， 当 预 判 到 遮挡 目标 时 ， 有 
这 时 目标 与 恒星 被 区 别 开 ， 且 只 对 目标 进行 跟踪 。 图 6d) Ne) 和 f) 的 前 一 帧 图 ， 图 6e) 和 
f) 是 对 同一 帧 图 的 不 同 跟踪 方法 ; 图 6e) 用 固定 波 门 跟踪 ， 波 门将 目标 和 恒星 同时 套 住 ， 图 
卡尔 曼 与 Camshif 结合 的 遮挡 目标 跟踪 方法 跟踪 目标 ， 这 时 只 对 目 


6f) 用 


变 光 目标 跟踪 
低 轨 遮挡 目标 跟踪 


图 6a) 和 b) 用 卡尔 曼 的 固定 波 门 跟踪 目标 ，a) 和 


偏离 波 门 中 心 ， 且 波 门 将 目标 和 恒星 同时 套 住 并 当 作 一 颗 更 大 的 上 
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b) 的 目标 分 别 在 波 门 的 右上 角 和 左 


标 跟踪 。 图 6c) 将 
日 本 文 的 遮挡 目标 跟踪 方法 跟踪 让 挡 目标 ， 


Fatt E 


标 跟 踪 时 的 遮挡 干扰 。 


与 Camshift 结合 的 遮挡 


6 ”遮挡 目标 跟踪 对 比 


标 跟踪 ， 消 除了 和 恒 


= 


图 7 的 横 坐 标 代表 目标 编号 ， 纵 坐标 代表 帧 数 ，z 的 位 置 代表 目标 编号 对 应 
了 直观 理解 ， 把 图 中 同一 种 方法 获得 的 各 目标 帧 数 连接 在 一 起 。 与 1 
目标 跟踪 捕获 的 遮挡 帧 数 及 捕获 遮挡 帧 


尔 曼 跟踪 相 比 ， 卡 尔 曼 
FP 的 有 效 帧 数 明 


显 高 于 卡尔 


曼 跟 踪 。 对 图 7 的 15 圈 目 标 进行 统计 ， 卡 尔 曼 跟踪 目标 时 谈 挡 帧 的 平均 捕获 率 为 45.7%， 


3.3.2 
图 


有 获 遮 挡 帧 中 平均 有 效 帧 率 为 30%; 卡尔 曼 与 Camshift 
获 讶 挡 帧 中 平均 有 效 帧 率 为 97.9%。 


低 罗 变 光 目 标 跟 踪 


时 ， 波 门 根据 目标 大 小 自 适应 变化 ， 收 敛 于 目标 质心 ,3 


Oo 


结合 的 遮挡 帧 3 


踪 效 果 。 


F 均 捕获 率 为 97.3%, 


对 应 的 图 8d). e) f) 


运动 速度 快 ， 目 标 


大 小 变化 频繁 的 目标 


的 6 帧 图 为 低 轨 快速 光 变 的 同一 目标 图 。 图 8a)、b)、c) 用 卡尔 曼 跟 踪 同 一 目标 ， 
当 目 标 大 小 频繁 变化 时 跟踪 波 门 固定 不 变 ，a) 由 于 目标 大 小 变化 频繁 
WIRE, B 8c) 目标 大 于 波 门 ， 不 易 准 确 捕 获 目 标 位 置 。 
Camshift 结合 的 自 适 应 的 可 变 波 门 跟踪 同一 目标 ， 当 遇 到 快速 运动 


| 卡尔 曼 与 


稳健 跟踪 上 


标 ， 最 终 获 


得 理想 的 跟 
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目标 编号 


注 : 深蓝 色 表示 每 圈 包 含 迹 挡 帧 数 ， 柳 色 表示 卡尔 曼 滤波 跟踪 捕获 帧 数 ， 灰 色 表 示 卡 尔 曼 滤波 跟踪 捕获 帧 中 的 有 
效 帧 ， 黄 色 表 示 本 文 改进 方法 捕获 帧 数 ， 浅 蓝 色 表示 本 文 改 进 方法 捕获 帧 中 的 有 效 帧 。 


7 ”遮挡 帧 的 跟踪 帧 数 和 有 效 帧 数 统计 


e) 
8 (RNA RE 


表 3 对 图 8 一 系列 无 拖 尾 的 低 轨 变 光 目 标 进行 捕获 率 统计 ， 由 表 可 知 卡尔 曼 与 Camshift 
结合 的 低 轨 变 光 空 间 目标 捕获 率 高 于 卡尔 曼 跟 踪 时 的 目标 捕获 率 ， 目 标 捕获 率 平均 提高 
18.4%。 


表 3 ” 低 轨 变 光 目 标 捕获 率 提升 
采集 时 间 HE FASA ” 波 门 大 小 是 否 TARUM MAK 


方法 . 和 $ Pa 

ae /min ” 轿 数 ”采集 帧 数 。 自 适应 变化 获 帧 数 /% 

卡尔 曼 3 15 207 ay 132 63.7 

卡尔 曼 与 Camshift 结合 3 15 207 是 170 82.1 


图 9 对 有 拖 尾 低 轨 短 曝 光 的 同一 空间 目标 进行 跟踪 ， 由 于 曝光 时 间 过 小 致使 除 目 标 外 
的 恒星 无 明显 拖 尾 。 图 9d). e f)s g)s hb). i) 是 对 a). bh c) 基于 概率 霍 夫 直线 检测 去 
除 目标 拖 尾 后 的 处 理 结 果 。 图 9d). eh f) 用 固定 波 门 跟踪 ， 当 目标 姿态 发 生变 化 时 波 门 保 
持 不 变 ， 由 于 目标 运动 速度 较 快 ， 图 9e) 帧 目标 被 跟 丢 ;与 图 9d), eh f) 对 应 的 g)、h)、 
i) 用 可 变 波 门 跟踪 变 光 目 标 ， 当 目标 姿态 发 生变 化 时 波 门 收敛 于 目标 。 

表 4 对 有 拖 尾 的 变 光 目标 进行 捕获 率 统 计 ， 卡 尔 曼 与 Camshift 结合 的 低 轨 变 光 空间 目 


> 
SI 
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| 


9 “有 拖 尾 的 低 轨 变 光 目 标 跟踪 


标 捕获 率 高 于 卡尔 受 跟 踪 时 的 目标 捕获 率 ， 目 标 捕 获 率 平 均 提 高 17.82%。 


表 4 有 拖 尾 的 变 光 目 标 捕获 率 提升 


方法 采集 时 间 ”目标 “平均 每 圈 ” 波 门 大 小 是 否 ”平均 每 圈 位 ”捕获 率 
/min 圈 数 ”采集 帧 数 。” 自 适应 变化 捕获 帧 数 /% 
卡尔 曼 3 10 230 T 138 60 
卡尔 曼 与 Camshift 结合 3 10 230 是 179 77.82 


图 10 对 图 8 和 图 9 的 一 类 目标 进行 跟踪 精度 统计 。 由 图 10 可 知 ， 卡 尔 曼 与 Camshift 
结合 的 跟踪 精度 优 于 卡尔 曼 跟 踊 时 的 跟踪 精度 ， 跟 踊 精 度 平均 提高 了 1.3 倍 ， 进 一 步 说 明 本 
文 提出 的 方法 的 可 行 性 。 
3.3.3 ”影响 目标 跟踪 的 因素 分 析 

针对 低 轨 变 光 目标 跟踪 过 程 中 经 常 出 现 薄 云 、 目 标 姿态 变化 等 导致 的 目标 星 像 丢 失 等 
青 况 ， 本 文 方法 不 具有 适用 性 。 由 于 不 能 检测 和 识别 目标 ， 所 以 更 不 能 跟踪 目标 ， 即 使 用 卡 
尔 曼 滤 波 的 预测 位 置 代替 Camshift 计算 出 的 目标 位 置 作为 观测 位 置 去 更 新 卡尔 曼 滤 波 ， 方 


= 
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123 4 5 6 7 8 9 10 1112 13 1415 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 


目标 编号 


法 
了 


起 作用 ;针对 


具有 一 定 的 适用 性 ， 但 仍 不 够 严谨 ， 因 
R 踪 数据 ， 数 据 误 差 也 相对 较 大 。 我 们 认为 : 针对 注 
， 然 后 进行 图 像 增强 ， 这 样 可 以 对 部 分 可 见 的 目标 进行 跟 


证: 蓝 色 表示 传统 卡尔 曼 ， 橙 色 表示 卡尔 曼 与 改进 的 Camshift 结合 。 


图 10 ” 变 光 目标 精度 提升 
为 预测 位 置 找 不 到 跟踪 对 象 来 跟踪 ， 
云 可 以 先 用 图 


pies 
LN » 


这 时 
像 二 值 化 和 滤波 去 噪 处 
但 对 星 像 丢失 的 目标 不 一 定 
图 像 增强 后 对 部 分 可 见 


UCET EN 


目标 ， 在 固定 时 


目标 姿态 变化 并 若隐若现 的 情况 ， 可 当 作 暗 弱 目标 进行 
的 目标 进行 跟踪 。 
捕获 率 是 对 某 一 圈 目 标 捕获 帧 数 与 该 轿 目 标 总 帧 数 的 比值 。 针 对 某 一 图 
间 曝 光 内 获得 的 目标 总 帧 数 是 不 变 的， 但 捕获 目标 位 置 的 帧 数 因 算 法 不 同 而 不 同 ， 当 捕获 
目标 位 置 的 帧 数 增加 时 ， 捕 获 率 将 提高 。 如 表 2、 表 3、 表 4 所 示 ， 与 卡尔 曼 跟 踪 相 比 ， 卡 


尔 曼 与 改进 的 Camshift 结合 的 每 圈 目 标 捕获 位 置 的 帧 数 有 所 增加 ， 捕 获 率 也 相应 增加 。 图 
像 处 理 时 间 指 卡尔 曼 与 Camshift 结合 
Camshift 跟踪 目标 所 用 的 时 间 做 对 比 ， 发 现 改进 的 Camshift 跟踪 目标 时 所 
即 程序 运行 效率 提升 。 图 


跟踪 目标 时 ，Camshift 跟踪 目标 所 


的 时 间 与 改进 的 


像 处 理 时 间 和 目标 捕获 率 存在 一 定 的 关系 : AR 


更 能 满足 
说 明 算 法 的 有 
当 目 标 与 


目标 跟踪 的 实时 性 需求 ， 使 得 目标 的 
图 像 处 理 时 间 无 太 大 关 


效 性 ， 与 


系 o 


恒星 星 像 


ESN, 


目标 与 4 


的 质心 发 生变 化 ， 这 时 计算 出 来 的 
目标 位 置 ， 也 属于 特例 ， 是 不 严谨 的 。 为 了 能 精准 和 稳定 跟踪 变 》 
目标 进行 跟踪 。 


部 分 的 有 效 


少 


将 卡尔 曼 与 Camshift 结合 ， 用 


4 结论 与 


空间 目标 与 碎片 观测 研究 对 于 实现 
意义 。 低 轨 变 光 


ary 
结合 :; 


Camshift 相 
稳健 跟踪 


目标 位 置 误差 过 大 ， 所 以 不 可 以 看 做 变 光 


h A> 时 


展望 


目标 跟 


颜色 提取 ， 降 低 颜 色 提 ! 


用 卡尔 曼 滤 波 外 推 
目标 ， 提 出 两 项 改进 ，(1) 将 Camshift 基于 彩色 图 


的 预测 值 代 蔡 观测 值 进行 星 像 遮 挡 的 


RAN Se 


EAE BL a AS KER) Ab, BE 
目标 ， 


时 间 变 短 了 ， 
像 处 理 时 间 变 快 时 
上 获 率 有 一 定 的 提升 ， 当 目标 的 捕获 率 提升 ， 


目标 
即使 捕获 
6 目标， 文中 


其 跟踪 识别 和 监测 预警 具有 重要 实践 价值 和 时 代 
踪 是 空间 目标 与 碎片 跟踪 的 重要 组 成 部 分 ， 本 文 将 卡尔 曼 滤 波 与 


预测 目标 运动 位 置 ，Camshift 跟踪 运动 


取 复 杂 性 ， 避 免 异 常 跟踪 ， 增 强 目 标 跟 踪 的 稳定 性 。 


标 预 判 方法 预 判 出 目标 与 恒星 是 否 存在 遮挡 ， 当 存在 遮挡 时 用 卡尔 曼 滤 波 的 预测 位 置 代替 


目标 。 为 确保 


像 的 颜色 提取 改 为 基于 gray 的 


(2) 用 遮挡 目 


4 期 王 恩 旺 ， 等 ， 基于 卡尔 曼 滤 波 与 Camshift 的 低 轨 空间 变 光 目 标 跟踪 方法 607 


Camshift 算法 计算 出 的 目标 位 置 ， 作 为 观测 值 去 更 新 卡尔 曼 滤波 。 

实验 结果 表明 ， 跟 踪 目 标的 波 门 由 固定 波 门 变 为 自 适应 的 可 变 波 门 ， 对 难 捕获 的 变 光 目 
标 具 有 很 好 的 跟踪 稳健 性 ， 对 瞬间 氮 挡 目标 获得 理想 跟踪 效果 ， 对 中 高 轨 目 标 跟踪 具有 普遍 
适用 性 ， 并且 增 强 了 跟踪 稳健 性 ， 提 高 了 精度 和 捕获 率 。 
本 文 方法 简单 ， 程 序 运 行 效率 快 ， 不 仅 适 用 于 单 目 标 稳定 跟踪 ， 也 适用 于 多 目标 稳定 跟 
踪 ， 而 且 自动 跟踪 目标 的 跟踪 效果 相对 理想 ， 有 科学 的 理论 支撑 ， 适 用 性 强 ， 还 能 为 国内 外 
空间 目标 跟踪 提出 一 种 方法 。 我 们 提出 的 方法 将 具有 广泛 应 用 前 景 ， 未 来 可 广泛 应 用 于 光电 
望远镜 阵 、 往 笛 望 远 镜 阵 、 光 学 巡天 望远镜 阵 和 大 口径 大 视 场 光学 望远镜 等 。 
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on Kalman Filter and Camshift 
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Abstract: In view of the frequent shape changes and fast moving speed of low-orbit dimming 
targets , it is difficult to accurately capture the target position information when tracking 
the target, a target tracking method combining Kalman filter and Camshift is proposed. 
According to the difficulty of capturing the target, Kalman filter extrapolation is introduced 
to predict the position of the target, tracking target with improved Camshift. That is, the 
target tracking method of single-channel white extraction based on gray image, to achieve 
the robust tracking of low-orbit variable light target and improve the target capture rate; 
When the target and the star are occluded each other instantaneously, the improved occluded 
target prediction method is used to predict the position. At the same time, the predicted 
position of the Kalman filter is used to replace the target position calculated by Camshift 
as the observation position to update the Kalman filter, so as to realize the robust tracking 
of the occluded target and improve the target acquisition rate. The experimental results 
show that the variable gate of the tracking target is adjusted adaptively according to the 
size of the target. It not only realizes stable tracking of medium-high orbit targets, but also 
has robust tracking effect for low-orbit dimming targets. This method has the advantages 
of fast running speed, high sensitivity, strong applicability and good real-time performance, 


and has good application value and broad application prospects. 


Key words: Kalman filtering; Camshift; dimming target; tracking 


